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ABSTRACT 

We use the direct Fourier meth od to calculate the redshift-space power spectrum 
of the maxBCG cluster catalog (jKoester et al.l [2007t ) - currently by far the largest 
existing galaxy cluster sample. The total number of clusters used in our analysis is 
12, 616. After accounting for the radial smearing effect caused by photometric redshift 
errors and also introducing a simple treatment for the nonlinear effects, we show 
that currently favored low matter density "concordance" ACDM cosmology provides 
a very good fit to the estimated power. Thanks to the large volume (~ 0.4 h^"^ Gpc'^), 
high clustering amplitude (linear effective bias parameter beS ~ 3 x (O.SS/trs)), and 
sufficiently high sampling density (~ 3 x 10~^ h"^ Mpc""^) the recovered power spectrum 
has high enough signal to noise to allow us to find weak evidence (~ 2a CL) for the 
baryonic acoustic oscillations (BAO). In case the clusters are additionally weighted by 
their richness the resulting power spectrum has slightly higher large-scale amplitude 
and smaller damping on small scales. As a result the confidence level for the BAO 
detection is somewhat increased: ~ 2.5a. The ability to detect BAO with relatively 
small number of clusters is encouraging in the light of several proposed large cluster 
surveys. 

Key virords: cosmology: observations - large-scale structure of Universe - galaxies: 
clusters: general - methods: statistical 



Since the flight of the Cobe satellite in the early 1990's 
the field of observational cosmology has witnessed ex- 
tremely rapid progress which has culminated in the estab- 
lishment of the St andard Model for cosmology - the "con - 
cordance" model ()Bahcall et al.lll999l : [Spergel et al.ll2003l ). 
This progress has been largely driven by the precise mea- 



Cosmic Microwave Background (CMB 


) jHananv et al.ll2000l: 


Netterfield et al.l |2002|: iBennett et all l2003l: iGrainge et al.1 


20031: Pearson et al.l 20031). However, i 


n order to break sev- 



eral degeneracies inherent in the CMB measurements one 
has to complement this data with other sources of informa- 
tion, such as the measurements of the SNe la lu minosity 
distances (|Riess et al.lll99 1 iPerlmutter et al.|[l999l ) or with 
the measurements of the large-scale structure (LSS) of the 
Universe as traced by galaxies or galaxy clusters. 

The simplest descriptor one can extract from the LSS 
measurements is the matter power spectrum. The broad- 
band shape of this spectrum is sensitive to the shape param- 
eter r = Qmh, and thus is useful in helping to establish the 
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currently favored low matter density "concordance" model, 
as well as the amount of baryons in relation to the total mat- 
ter fb — Qb/^m- Currently the two largest galaxy redshift 
surveys are the 2dF Galaxy Redshift SurvejQ (2dFGRS) 
and the Sloan Digital Sky Survej0 (SDSS) with its latest 
data release 5 (DR5), providing redshifts to ~ 220,000 and 
~ 675, 000 galaxies, respectively. The SDSS galaxy sample 
consists of two broad classes: (i) the MAIN galaxy sample, 
reaching redshifts of z ~ 0.25; (ii) the Luminous Red Galaxy 
(LRG) sample covering redshift range z ~ 0.15 — 0.5. Some 
other characteristics of these surveys are listed in Table [T] 

In addition to the precise measurement of the broad- 
band shape of the power spectrum, SDSS LRG and 2dFGRS 
samples have proven useful in allowing the detection of theo- 
retically predicted oscillatory features in the spectrum. The 
source of these spectral fluctuations is the same as that giv- 
ing rise to the prominent peak structure in the CMB angular 
power spectrum - acoustic oscillations in the tightly cou- 
pled baryon-photon fluid p rior to the epoch o f recombina- 
tion (jSunvaev fc Zeldovichiil970l : [Peebles fc Yu|[r97Q ). How- 



^ http:/ /www. mso.anu.edu.au/2dFGRS/ 
^ http://www.sdss.org/ 
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Table 1. Some characteristics of the 2dFGRS, SDSS DR5, and maxBCG samples. 



Survey 


Number of 


Sky area 


Redshift 


Comoving volume 


Effective bias parameter wrt to 




objects 


{deg2) 


coverage 


(h-^Gpc^*) 


the model with erg = 0.85 


2dFGRS 


~ 220, 000 


~ 1200 


£ 0.3 


~ 0.07 


~ 1.0 


SDSS DR5 


~ 675, 000 


~ 5700 


MAIN: < 0.25 


MAIN: ~ 0.2 


MAIN: ~ 1.2 








LRG: ~'^.15 - 0.5 


LRG: ~ 1.3 


LRG: ~ 2.0 



maxBCG ~ 13, 800 ~ 7000 0.1 - 0.3 ~ 0.4 ~ 3 



ever, the relative level of fluctuations in the matter power 
spectrum is strongly reduced compared to the correspond- 
ing fluctuations in the CMB angular spectrum, owing to the 
fact that most of the matter in the Universe is provided by 
the cold dark matter (CDM) that does not participate in 
acoustic oscillations. After recombination, or more precisely 
after the end of the so-called "drag-epoch" which happens 
at redshifts of a few hundred, as the baryons are released 
from the supporting pressure of the photon gas they start 
to fall back to the CDM potential wells that have started to 
grow already since the matter-radiation equality. Since the 
amount of baryons in the total matter budget is not com- 
pletely negligible - baryons make up roughly one fifth - the 
acoustic structure imprinted onto the baryonic component 
also gets transferred to the spatial distribution of the CDM. 
After baryonic and CDM components have relaxed one ends 
up with ~ 5% relative fluctuations in the total mater power 
spectrum. 

In order to observe these relatively small fluctuations 
one needs to have: 

(i) Large survey volumes, to reduce cosmic variance; 

(ii) High enough spatial sampling density, to decrease dis- 
creteness noise. 

Evidently, for the fixed telescope time there is a trade-off 
between these two requirements. Amongst the available 
spectroscopic samples the optimum is currently best met 
by the SDSS LRGs, and indeed, at the beginning of 2005 
the detection of the acoustic peak in the spatial two-point 
correl ation function was announced by lEisenstein et al.l 
l|2005h . At the same time the final power spectrum mea- 
surements of the 2dF GRS also reveal ed the existence of 
the acoustic features (|Cole et al.ll2005l ). Baryonic acoustic 
oscillations (BAO) in the power spectrum of the SDSS 
LRG (DR4) sample were found by iHiitsil l|2006al fl and 
the correspondi ng cosmologic al parameter estimation was 
carried out in iHiitsil (|2006d '). Also, several more recent 
detect ions of BAO have been reported: iPadmanabhan et al.l 



l|2007l ) a ndlBlake et"al] ll2007f) u se d SDSS photomet ri c LRG 



sample, iTegmark et al 



Martfnez et al.l 



jgOOef ) : lOkumura et al.] 



magal (u 



Cabre fc Gaztafiagal ll2009f): iGaztanaga et al 



200S 



(200g 



Sanchez et all (|2009l ): iKazin et al 



20091) analyzed spectroscopic LRG sample, iPercival et al 
2007g) carried out a combined analysis of the SDSS 



MAIN and LRG spec t roscop ic galaxy samples, and finally, 
IPercival et all l|20[)7al . |2009| ) analyzed SDS S LRG, MAIN 
and 2 dFGRS samples. In the recent paper bv lNomura et all 
l|2009l) a weak detection of the BAO damping is reported. 
The usefulness of BAO arises from the fact that it can 

^ For an independent comparison see lYamamoto et all fcOOSl ) . 



provide us with a standard ruler in the form of the size of 
the sound horizon at decoupling, enabling us to carry out a 
purely geometric cosmological test: comparing the apparent 
size of the ruler along and perpendicular to the line of sight 
with the physical size of the sound horizon (which can be 
well calibrated from the CMB data) one is able to find Hub- 
ble parameter H{z) and angular dia meter dist ance dAiz) 
corresponding to the redshift z (e.g. IHu fc HaTn ian 2003). 
Obviously, having determined H{z) and dA{z) one can put 
constraints on the dark energy (DE) equation of state pa- 
rameter wde{z). In that respect radial modes, which enable 
us to find H{z), are more useful since H{z) is related to 
wde{z) through a single integration whereas dA(z), which 
is given by the angular modes, involves double integration 
over wde{z). Thus one would certainly benefit from accu- 
rate redshift information. 

If one wishes to use the BAO signal as a precise standard 
ruler there are a few complications one has to overcome: 

(i) As the density field goes nonlinear couplings between 
Fourier modes (and thus the resulting transfer of power) 
will modify the BAO signal, leading to the damping of the 
oscillations; 

(ii) In order to fully exploit the information in 
galaxy /cluster surveys one has to understand the relation 
of these objects to the underlying matter density field, i.e. 
one has to have a realistic model for biasing. 

These two complications can be fully addressed only 
through costly N-body simulations. There are several pa- 
pers that have investigated the detectability and possible 
systematics of BAO extrac t ion u si ng the numeri c al sim - 
ulations e.g. iMeiksin et al.l (Il999t): ISpringel et all ll2005h: 



ISeo fc EisensteinI ||2005|): iHuff et al.1 (12 007'): 'Koehl er et all 
(|2007l ): lAngulo et alj (|2008l '): ISmithet al. (2008). However, 
significantly more work towards these directions is probably 
needed before we can harness the full power provided by the 
potentially very clean geometric tool in the form of BAO in 
the matter power spectrum. 

It is worth pointing out that actually one of the first 
possible hints for the existence of BAO in the matter power 
spectrum came from the a nalysis of the spat ial clustering 
of the Abell/ACO cluster s (iMiller et al.ll200lM ^ The power 
spectrum measurement in lMiller et al.l I 200lf ) used 637 clus- 
ters with richness R ^ 1. In this paper we analyze the 



^ To be more precise again, at the end of the "drag-epoch". In 
the "concordance" ACDM model the sound horizon at the "drag- 
epoch" is ~ 5% larger compared to the one at recombination. 
® Interesting features in the spatial two-point correlation func- 
tion of the Abell/A CO cluster sample were also revealed by 
lEinasto et"all lll997l '). 
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Figure 1. The effective volume (see Eq. ©) of the maxBCG Figure 2. Angular distribution of 12,616 maxBCG clusters in 

cluster sample in comparison to the 2dFGRS and SDSS MAIN SDSS survey coordinates {rj, A). The 5x5 chessboard pattern of 

and LRG samples from data release five (DR5). black/gray points shows the angular part of the division used for 

"jackknife" error analysis. The union of solid rectangles represents 
our reconstruction of the angular mask. 



maxBCG cluster sample (jKoester et al.l[2007h which is cur- 
rently by far the largest cluster catalog available, spanning a 
redshift range z = 0.1 — 0.3, containing 13, 823 objects, and 
covermg ~ OAh-^ Gpc^ of comoving volume (see Table [T] 
for comparison with 2dFGRS and SDSS DR5). The biggest 
advantage of using galaxy clusters instead of the typical, i.e. 
L*-galajcies, is the fact that their spatial clustering signal 
is strongly amplified with respect to the clustering of the 
underlying matter: bias parameter b often reaching values 
~ 3 and above. In Fig. [T] we have compared the effective 
volume (see Eq. ([2])) of the maxBCG survey with the two 
other surveys that have lead to the detection of BAO: 2dF- 
GRS and SDSS LRG. In addition we have also shown the 
effective volume of the SDSS MAIN galaxy sample. Under 
the assumption of Gaussianity the power spectrum measure- 
ment errors AP are simply found as: AP/P = 2/Vefs/Vk, 
where 14 is the volume of the shell in fc-space over which 
the angular average is taken, i.e. Vk = 47rfc^Afc/(27r)^. It 
is no surprise that the SDSS LRG sample is currently un- 
beatable. However, we notice that maxBCG sample should 
provide better measurement of the power spectrum on scales 
larger than k ~ 0.1 /iMpc"^ compared to the 2dFGRS. The 
rapid fall of the effective volume of the maxBCG sample on 
smaller scales is caused by t he photometric reds hift errors 
of Sz ~ 0.01 as estimated bv lKoester et all {200t) . 

Thus one would expect the maxBCG sample to reveal 
acoustic features. As our further analysis shows this indeed 
turns out to be the case, albeit with a relatively modest 
confidence level. 

Our paper is organized as follows. In Section 2 we give 
a brief description of the maxBCG catalog with the corre- 
sponding selection effects. The core part of this paper, which 
is devoted to the power spectrum analysis, is given in Sec- 
tion 3. In Section 4 we test the stability of our spectrum 
measurements by relaxing several assumptions made in the 
original analysis. Finally, Section 5 contains our conclusions. 



2 DATA AND SURVEY SELECTION 
FUNCTION 

A maxBCG cluster catalog l|Koester et al.l |2007| ) is com- 
piled via the "r ed sequence" cluster detecti on method 
'^Ostrander et al.| [l998: Gladdcrs & Ycc 2000; B ahcall et all 
120031 ) applied to the SDSS photometric data. The catalog 
contains 13, 823 galaxy clusters with velocity dispersions 
greater than ~ 400 km/s and covers ~ 7000 square de- 
grees of sky between redshifts 0.1 and 0.3. The photometric 
redshifts are estimated using the tight relation between the 
ridgeline color and redshift and are claimed to have accuracy 
Sz = \/ ((zphoto — Zapec)^) ~ 0.01 esseutially independent of 
redshift. For more details about the catalog we refer to the 
original source [Koester et al.l (|2007l ). 

In our power spectrum analysis we are going to neglect 
the three southern SDSS stripes, leaving us with 12, 616 
galaxy clusters over ~ 6800 square degrees of sky. With 
the redshift range 0.1 < z < 0.3 this corresponds to 
~ 0.4 Gpc'^ of comoving cosmic volume. Neglecting the 
three narrow southern stripes helps us to achieve a better 
behaved survey window function with reduced sidelobes, and 
thus with reduced leakage of power. 

The angular distribution of the remaining clusters is 
given in Fig. [2] Here the angular coordinat es are plotted us- 
ing th e SDSS survey coordinates (77, A) (e.g. IStoughton et al.l 
l2002f ). The chessboard pattern of black/gray points in Fig. 
[2] represents the division of the sky into 5x5 angular re- 
gions used for "jackknife" error analysis. We also use three 
divisions along the redshift direction, making the total of 

3 X 5 X 5 = 75 regions, each containing approximately 168 
clusters. The union of solid rectangles in Fig. [2] represents 
our reconstruction of the survey angular mask. All the rect- 
angles are aligned with the SDSS imaging scan stripes. As 



^ http:/ /umsdss. physics. lsa.umich.edu/catalogs/ 
maxbcg_public_catalog.dat 



4 Gert Hiitsi 



900 



o 

CO 

?3 800 - 




100 I 1 ' 1 1 1 1 ' 1 1 1 

0.1 0.12 0.14 0.16 0.18 0.2 0.22 0.24 0.26 0.28 0.3 

redshift z 

Figure 3. Redshift histogram of the maxBCG sample. The solid 
smooth line corresponds to the cubic spline fit and the dashed 
line represents the expected number of clusters per redshift bin 
for the volume- limited sample. 

the number density of galaxy clusters in the maxBCG sam- 
ple is rather high, one can determine relatively accurately 
the beginning, ending, and also possible gaps in the scan 
stripes. We call the angular mask obtained this way a "min- 
imal" mask in contrast to the "maximal" mask, which is 
built in the same manner except that each of the rectangles 
is extended by an amount corresponding to the mean cluster 
separation. We carry out our power spectrum analysis using 
both of these angular masks. 

In Fig. Owe show the redshift distribution of maxBCG 
clusters. Here the solid smooth line corresponds to the cubic 
spline fit and the dashed line represents the expected num- 
ber of clusters per redshift bin Az — 0.2/30 for the volume- 
limited samplqj assuming a spatially flat ACDM cosmology 
with flrn = 0.27. To convert this redshift distribution to 
the radial selection function n(r), i.e. the comoving number 
density of clusters at distance r, we fix the background cos- 
mology the low matter density model mentioned above. Our 
power spectrum calculations are done for both of the radial 
selection models shown in Fig. [S] 

For the full survey selection function n(r) we assume, 
as usual, that it can be factorized to the product of angular 
and radial parts, i.e. n(r) = n{r)n{r) (r = |r|, f = r/r). 
Here n(r) is the radial selection function and the angular 
selection part n(f) is assumed to take value 1 if the unit 
direction vector f is inside the survey mask and otherwise, 
i.e. we assume a simple binary angular mask. 

3 POWER SPECTRUM ANALYSIS 

We calculate the redshift-space power spectrum of the 
maxBCG catalo g using the dir ect Fourier method as de- 
scribed in iFeldman et all (|l994i l (FKP). Strictly speaking, 

^ In iKoester et al.l it is suggested that maxBCG sample 

is rather close to being volume-limited. 



power spectra determined in this way are so-called pseu- 
dospectra, meaning that the estimates derived are convolved 
with a survey window. Since in the case of the analyzed 
maxBCG sample the volume covered is very large, reach- 
ing OAh~^ Gpc'^, and also the survey volume has relatively 
large dimensions along all perpendicular directions, the cor- 
relations in the Fourier space caused by the survey window 
are rather compact. On intermediate scales and in the case 
the power spectrum binning is chosen wide enough, FKP es- 
timator gives a good approximation to the true underlying 
power. 

Instead of dir ect summation as presented in 
iFeldman etahl l| 19941 ) we speed up the calculations using 
Fast Fourier Transforms (FFTs). This gives rise to some ex- 
tra complications. As our density fiel d, which is built using 
the Triangular Shaped Cloud (TSC) l|Hocknev fc EastwoodI 
1988) mass assignment scheme, is now "restricted to live" 
on a regular grid with a finite cell size, we have to correct for 
the smoothing effect this has caused. Also, if our underlying 
density field contains spatial modes with higher frequency 
than our grid's Nyquist frequency, /cNy, then these will 
be "folded back" into the frequency interval the grid can 
support, increasing power close to fcNy - the so-called 
aliasing effect. Often there is no need to correct for the 
aliasing effect since one usually chooses fine enough grid 
that has Nyquist frequency significantly higher than the 
highest spatial frequencies of interest. Although the aliasing 
has a very mild effect on our results w e cor r ect it by using 
the iterative scheme as described in with a 

slight modification: we do not approximate the small-scale 
spectrum by a simple power law, but also allow for the 
possible running of the spectral index, i.e. the parametric 
shape of the power spectrum is taken to be a parabola in 
log-log. We have co mpare d our antialiasing scheme with 
the original recipe of Ijind |2p05) finding a full agreement 
between the two. Thus the measured spectrum along with 
the acoustic features is very insensitive to the variations in 
the antialiasing method applied. 

For the full details of our power spectrum calcu lation 
method along with several tests we refer the reader to Hi itsil 
(|2006bl lal). As a very brief summary, our spectrum determi- 
nation consists of the following steps: 

• Determination of the survey selection function, i.e. 
mean expected number density without any clustering n(r) 
(see Section [2}. This smooth field, with respect to which 
the fiuctuations are measured, is modeled using a random 
catalog that contains 100 times more objects than the real 
maxBCG catalog, i.e. 1, 261, 600 in total; 

• Calculation of the overdensity field on a 512^-grid using 
TSC mass assignment scheme; 

• FFT of the gridded overdensity field; 

• Calculation of the raw 3D power spectrum by taking 
the modulus squared from the output of the previous step; 

• Subtraction of the shot noise component from the raw 
3D spectrum; 

• Recovery of the angle averaged spectru m usin g a mod- 
ified version of the iterative scheme of iJinj l|2005h . 

The results of our power spectrum calculation along 
with fitted models are shown in Figs. |4] and [5] Here the 
errorbars AP are calculated as described in FKP, which 
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Figure 4. Upper panel: Filled circles with errorbars show the 
redshift-space power spectrum of maxBCG cluster sample. Solid 
line represents the best fitting model in case of three free pa- 
rameters: bgff, q, a, and assumes Monte Carlo error model. Long 
dashed line shows the same model without survey window con- 
volution applied and short dashed line corresponds to the linear 
model spectrum. Lower panel: the same as above, with the broad- 
band shape of the spectrum removed by dividing with a "smooth" 
model spectrum without BAO. 



assumes that the density fie ld follows Gaussian statistics 
(see also lTegmark et al.|[l99l ): 



AP _ 



where the effective volume VcEf (see Fig. [T| is given by: 



I 



n{r)P{k,ii) 



(1) 



(2) 



and Vfc is the volume of the fc-space shell over which the 
angular average of the 3D spectrum is taken, i.e. 



Vk 



(27r)3 



(3) 



In Eq. rniin and rmax are comoving distances correspond- 
ing to the redshift bounds of the survey (0.1 and 0.3 in case of 
maxBCG), n{r) is the radial selection function, and Q is the 
area of the sky in steradians covered by the survey. P{k, is 
the anisotropic power spectrun{f| where ^ is the cosine of the 
angle between the line of sight direction r and wavevector 



Here we assume a flat sky approximation. 
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Figure 5. Same as Fig. |4l with the only difference that here we 
have two free parameters: b^g, q, and the value for a is fixed to 

ao/i-^Mpc. 



k, i.e. fj, = r li. To calculate effective volumes for the spec- 
troscopic surveys, such as 2dFGRS and SDSS shown in Fig. 
[1] one can neglect mild distortions due to redshift-space dis- 
tortions and take power spectrum to be a function of k only, 
i.e. drop integral over /i in Eq. ((2)1. However, for photometric 
surveys with large distortions along the line of sight this ap- 
proximation is no longer valid. For the maxBCG survey we 
fit a parametric model in the form of Eqs. (|6]) and ((T} to the 
measured isotropized power spectrum. To calculate effective 
volume plotted in Fig. [T]we have taken P(fc,/x) as in Eq. 
(O with the replacement f{k) — > f{k,ii) = exp[—{fika)^]. 
Here a is the spatial smoothing scale along the line of sight; 
a — S-Sz . 



the power spectru m variance via the "jackknife" method (see 
e.g. lLuptonlll993l ). For this purpose we have divided the sur- 



In addition to the FKP errors we have also estimated 
; power spectru r 
. lLuptonlll993l l. 
vey volume into 5x5 angular regions (see Fig. [2]) along with 
3 redshift bins, making the total of 75 regions, each contain- 
ing on average 168 clusters. Now the power spectrum is cal- 
culated 75 times with each time one of the regions omitted. 
The covariance of the power spectrum can now be estimated 
via the scatter of the power spectrum measurements: 



COv(Pi,Pj) = 

Pi = 

Here Pj'"' represents the power spectrum measurement for 
the i-th power spectrum bin in case the n-th survey subvol- 



1 ^ 



iV 



(5) 
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Figure 6. Example slices (5 degrees wide) through the mock 
catalogs and maxBCG sample (bottom right-hand panel). All of 
the point sets have very similar two-point functions. 



ume was excluded from the calculations. There are several 
obvious problems related to the "jackknife" procedure, e.g. 
how to choose the size for the subvolumes, also the num- 
ber of subvolumes in our case is relatively small, leading to 
a rather noisy covariance matrix. To overcome those short- 
comings we have used a third method to estimate the er- 
rors. In this Monte Carlo approach we have generated 1000 
mock catalogs with the survey selection function matching 
that of the maxBCG sa mple u s ing a m ethod described in 
detail in Appendix B of iHtitsil (|2006al ) . Our method con- 
sists of three parts: (i) generation of the density field using 
the optimized 2nd order Lagrangian perturbation scheme 
(2LPT), (ii) Poisson sampling of the density field with the 
sampler whose parameters are tuned to result in the clus- 
tering strength similar to the maxBCG sample, (iii) addi- 
tion of the redshift-space distortions and photo-z errors. 
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Figure 7. Comparison of the FKP, "jackknife", and Monte Carlo 
error estimates. Note that the errors as given here were calculated 
by dividing each mock spectrum by a smooth fit to it, i.e. the 
additional uncertainty due to the variation of the power spectrum 
amplitude is taken out. 

In Fig. [6] we show a few five degrees wide slices through 
mock catalogs along with the slice through the maxBCG 
sample (bottom right-hand panel). All of the point sam- 
ples in Fig. |6] have very similar two-point functions and 
since the 2LPT is able to capture the morphology of the 
quasilinear density field down to scales of ~ Wh~^ Mpc we 
expect also a rather good match for the higher order cor- 
relations on those scales. All of the three error estimates 
are compared in Fig. [7] where for clarity we have removed 
the smooth broad-band component of the spectrum. We see 
that in general errors are in good agreement. In particu- 
lar, the Monte Carlo and FKP errors agree remarkably well 
on almost all, except for the very largest scales. The mild 
discrepancy at the largest scales is probably due to the fi- 
nite box size (which leads to a poor sampling of the very 
largest Fourier modes) used in our 2LPT simulations. In 
Fig. [S] we compare the "jackknife" and Monte Carlo corre- 
lation matrices corr(Pi,Pj) = cov(Pi, Pj)/^var(Pi)var(Pj). 
Due to relatively small number of subvolumes (75) used in 
our "jackknife" analysis the resulting correlation matrix is 
quite noisy. Also, one can notice that on the largest scales 
the matrix elements neighboring the main diagonal drop off 
quite a bit faster in comparison to the Monte Carlo case. 
In the following model fitting part of our paper we use all 
three error estimates obtained above. For the case of the 
FKP errors a simple diagonal covariance matrix is assumed. 
If not stated explicitly, all of the figures in this paper assume 
Monte Carlo covariance matrix. 

Before we can fit model spectra to the data there are 
several effects one has to take into account: 

• Photometric redshift errors lead to significant damping 
of the spectrum; 

• Our observed spectrum is a pseudospectrum and thus 
before fitting one has to convolve model spectra with the 
survey window; 

• On smaller scales nonlinear effects become noticeable 
and one needs a model to account for these; 
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Figure 8. "Jackknife" (upper panel) and Monte 
Carlo (lower panel) correlation matrices corr(Pi,Pj) = 
cov(Pi,Pj)/^var(Pi)var(Pj). 



• And finally, even without ptiotometric redsliift errors 
there are spectral modifications caused by the redshift-space 
distortions. 

Our model for the relation between the linear spectrum 
Piin(fc) and the observed one Pobs(fc) that tries to accom- 
modate all the four effects listed above reads as: 



Pobs(fc) = / dk' k'^P{k')K{k',k) , 



where 



J'(fc) = &cff (i + gfc^)/(fc)Piin(fc). 



(6) 



(7) 



Here boff is the cluster bias parameter that also incorpo- 
rates the boost of the power at large scales due to the linear 
redshift-space distortions. Smaller scale redshift distortions 
are fully eliminated if the cluster detection algorithm is able 
to perform a perfect finger-of-God compression. However, in 
reality the compression is not ideal and thus there should be 
some extra uncertainties for the central redshift of the galaxy 
cluster. This can be modelled as an additional smoothing ef- 
fect in addition to the photometric redshift errors of the indi- 
vidual galaxies. The factor (l + qk^^^) is our simplistic model 



to treat effects due to nonlinear evolution. This fo rm is very 
si milar to the o ne suggested in ICole et all l|2005h . however 
m iHiitsil (|2006d ) we found that thepower law index 3/2 pro- 
vides better fit to the Halo Modejj results compared to the 
value 2 used in lCole et all l|2005l ). The function f{k) in Eq. 
(O models the damping of the spectrum due to photometric 
redshift (photo-z) errors. Under the flat sky approximation 
it is easy to derive the analytic form for f{k) which reads: 



/(.) = ^erf(.k)^ 



(8) 



where the error function erf(x) = /^^ exp(— t^) dt . This 
result assumes that photo-z errors follow Gaussian distribu- 
tion with dispersion Sz and a is the corresponding spatial 
smoothing scale, i.e. a = -^Sz . And finally, Eq. ^ models 
the convolving effect of the survey window. We find that the 
mode coupling kernels K{k' , k) can be well described by the 
following analytic fit: 



K{k',k) = K{k,k') 
where 



^[gik + k')^g{k-k')] 



g{x) — arctan 



(- 



(9) 



(10) 



and the normalization constant C is derived by demanding 
that / K{k, k')k''^ dk' = 1, giving: 



C : 



-iih\ 2 



(11) 



In the case of the maxBCG survey geometry the best fit- 
ting parameters a and b axe found to be: a = 0.0044 and 
b — 0.0041. For the motivation of this analytic form and 
for the complete description on how these kernels are de- 
termined see [Hiitsi (2006^) . The analytic approximation is 
compared with the numerically determined coupling kernels 
in Fig.[9l We see that our analytic form provides a very good 
fit to the core parts of the numerical kernels with only slight 
deviations at the wing portions. In practice, performing the 
convolution of the model spectra with those approximate 
kernels in place of the precise numerical ones leads to com- 
pletely negligible differences. 

In this paper we keep the background cosmology fixed 
to the "concordance" flat ACDM model with Q,rn = 0.27, 
Slj, = 0.045, and h = 0.7. The bias parameters quoted below 
are measured with respect to the model with erg — 0.85. Our 
decision not to fit for the cosmological parameters is based 
on two main reasons: 



(i) Photo-z uncertainties that are given in iKoester et al.l 

are tested only for the case of the brightest cluster 
members that have spectroscopic redshifts available. The 
errors for the whole galaxy population in clusters might 
differ, and be possibly larger than the numbers given in 
iKoester et al] (2007); 

(ii) Uncertainties in the cluster finding algorithm which 
are hard to quantify and which might lead to the additional 



^ For a comprehensive review on Halo Model see lCoorav &: ShethI 
1I2OO2I) . 

A simpler analytic form with a = b = 0.0041 gives almost as 
fit. 



8 Gert Hiitsi 




0.15 0.2 
k [hMpc^i] 

Figure 9. Mode coupling kernels. Dashed lines show the analytic 
approximation given in Eqs. JToJ-JTlJ. 
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Figure 10. Error ellipses for the case with three free param- 
eters. The solid, dashed, and dotted lines correspond to the 
Monte Carlo, "jackknife" and FKP error models, respectively. 
The crosses mark the best fitting parameter values. Note the 
strong degeneracy between fe^g and a. 



smoothing/overmerging, and thus to the extra uncertainties 
in derived cluster redshifts. 

Any untreated systematics in the determination of the ra- 
dial smoothing scale a will start to interfere with the cosmo- 
logical parameter estimation, resulting in biased estimates. 
Due to the uncertainties listed above we first treat the ra- 
dial smoothing scale a as a free parameter. Thus we are left 
with three free parameters: 6eff, o", q. We also investigate 
separately the case where the radial smoothin g scale is fixed 
at c/H q ■ 5z ~ 30/i"^Mpc as suggested in iKoester et al.l 
Some of the results of these calculations are pre- 
sented in Figs. |4] and O and in Tables [2] and [S] To calculate 
theoretical model spectra we have used fitting formulae fo r 
the transfer functions presented in lEisenstein fc Hul l|l998l ). 
In the lower panels of Figs. |4] and [5] we have divided the 



spectra by the "smooth " model spectrum with acoustic os- 
cillations removed (see lEisenstein fc Hul [l99^ for details). 
For fitting purposes we have used Levenberg-Marquardt x^- 
minimization technique as described in IPress et al.l (Il992l ) 
with the additional modification to allow for the nondiag- 
onal data covariance matrices. We perform the fitting for 
all three error models: (i) FKP errors with diagonal covari- 
ance matrix, (ii) "jackknife" and (iii) Monte Carlo errors 
along with full covariance matrices. The solid lines in Figs, 
m and [5] represent the best fitting "wiggly" models for the 
case of Monte Carlo errors. In these figures we also demon- 
strate the damping effect of the window convolution. The 
corresponding linear spectra are shown with short dashed 
lines. One can see that the shape of the measured spectrum 
deviates strongly from the linear model spectrum, which is 
mainly driven by the photo-z errors. The best fitting pa- 
rameter values along with the corresponding values are 
given in Tables [21 E] There we also present the x^ value for 
the best fitting "smooth" model. In general we can see that 
the results for the x^-statistic agree well with expectations 
13.0±5.1 and 14.0±5.3 for 13 and 14 degrees of freedom, rele- 
vant for Table[2]and|3l respectively. Comparing the x^ values 
for the best fitting "smooth" and "wiggly" models one can 
assess the confidence level for the BAO detection, the results 
of which are presented in the last columns of Table [2] and |3| 
We note that using the full "jackknife" and and Monte Carlo 
covariance matrices in place of a simple FKP error model 
we obtain higher confidence for the detection of the oscilla- 
tory features in the spectrum. This is easy to understand: 
our first error model, where we assume diagonal covariance 
matrix and FKP errors is certainly not realistic. Survey win- 
dow, and more importantly nonlinear mode-mode coupling 
induces significant correlations between neighboring power 
spectrum bins. Stronger correlations amongst the bins make 
the spectrum "harder to distort" , and thus any surviving os- 
cillatory features gain higher statistical weight. It is interest- 
ing to note that in case we allow cr to be a free parameter the 
best fitting values are certainly larger than 30 Mpc. Also 
the bias parameter feeff has a noticeable increase compared 
to the cr = 30 Mpc case. In fact there is a strong degen- 
eracy between b^s and a: if one chooses to increase 6ofi, one 
can still obtain a good fit to the data by also increasing a by 
a suitable amount (see Fig. llOp . The fact that the case with 
three free parameters favors higher value for a (and thus also 
for fooff ) is mainly driven by the relatively large amplitudes of 
the two first power spectrum bins. On the other hand, if one 
assumes that a = 30 Mpc is a reasonably correct value 
then one is left with some discrepancy on the largest scales 
between the model and data (see Fig. [5}. It is interesting 
that some large-scale "extra p ower" has also been detecte d 
by several other a uthors e.g. iPadmanabhan et al. I l|2007l ): 
iBlake etall (|20o3) who analyzed SDSS photometric LRG 
sample^fof course these mild discrepancies on the largest 
scales might be due to inaccurate treatment of the survey 
selection effects. Taking into account several uncertainties 
concerning the appropriate value for a, and also uncertain- 
ties concerning the survey selection effects, we are here un- 
fortunately unable to draw any firm conclusions about the 



" See also Fig. 4 in lHirtsill2006al where SDSS spectroscopic LRG 
sample reveals some mild excess of power on the largest scales. 
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Table 2. Best fitting model parameters and values along with the inferred confidence levels for BAO detection for the case with three 
free parameters: feeffi ?! f''- 
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BAO detection 
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Table 3. Analog of Tabic [2] for the case with two free parameters: b^ff i 9i a^nd the value for cr is fixed to 30/i ^ Mpc. 
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"Jackknife" 


5.70 ±0.37 


12.1 ±2.4 


12.4 


15.7 
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Monte Carlo 
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existence of the excess large-scale power. This issue certainly 
deserves a dedicated study and currently the best available 
sample for this purpose is probably the SDSS LRG spectro- 
scopic sample. 

However, one extra exercise we can perform here is to 
try to investigate how well the obtained values for the bias 
parameter beff agree with the model expectations. If one as- 
sumes a flat ACDM model with the parameters as given 
above (i.e. ilm = 0.27, as = 0.85) then for the mass-limited 
cluster survey covering the same redshift range and sky area 
as maxBCG sample, and giving the same number of clusters, 
one infers a lower mass limit of Miow = 7.1 x lO" H'^Mq. 
We h ave calculated Mipv, usin g the Sheth-Tormen mass func- 
tion (Shcth & Tormen"l999l) averaged over the past light- 
cone (sec Hiitsi 2006b for details). The corresponding light- 
cone averaged effective bias parameter including increase 
of po wer due to la rge-scale redshift-space distortions (see 
again iHiitsil I2OO6IJ ) with respect to the z = linear spec- 
trum turns out to be feoff — 2.8 (b^g ~ 7.8)13 Comparing 
this value of ~ 8 with the ones given in Tables [2] and 
[3] one sees that a = 30/i~^Mpc case has somewhat lower 
bias value, whereas the case with free a tends to have val- 
ues for 6eff, which seem to agree reasonably well (especially 
true for the case with "jackknife" errors). Thus this simple 
analysis would favor the "extra smoothing" hypothesis over 
the "extra power" one. 



In reality of course for the mass-limited sample the number 
density drops with redshift, while maxBCG catalog is close to 
being volume-limited. We have calculated that to achieve a co- 
moving number density of 3.2xl0~^ h'^ Mpc"'', as relevant for the 
maxBCG catalog, the lower mass limits for redshifts 0.1 and 0.3 
should be 7.8 X lO^^ H'^Mq and 6.6 X lO^^ H-^Mq, respectively. 
The effective bias parameters with respect to the z = linear 
matter power spectrum for these redshifts and lower mass limits 
are 6Jff°'^ — 2.75 and b^g^'^ — 2.78, i.e. basically independent of 
redshift, and also agreeing very well with the light-cone averaged 
value of 2.8 quoted in the main text. This approximate constancy 
of the clustering amplitude at different redshifts for the objects 
with fixed comoving number density is a well known result. 



4 STABILITY OF THE RESULTS 

In this Section we are going to relax several assumptions 
that were made in our original analysis in order to test for 
the stability of the BAO features. The results of these inves- 
tigations are presented in Fig. [11] We start by testing the 
stability of the BAO against possible uncertainties in the 
survey selection function. For this purpose we have repeated 
power spectrum calculations with modified radial and angu- 
lar selections. On the upper left-hand panel of Fig. [11] we 
show the change induced in the BAO if one replaces the 
smooth cubic spline fltted radial selection with the one cor- 
responding to the volume- limited sample (see also Fig. [3]). 
Here the solid line represents the case with a new radial se- 
lection whereas the dashed line with errorbars corresponds 
to our original BAO measurementQ We see that at almost 
all the scales, except for the very largest ones, the BAO fea- 
tures are essentially unchanged. However, this exercise shows 
that the largest scales are very vulnerable with respect to 
the unaccounted slow trends in the survey selection. The 
results of the similar exercise for the angular selection func- 
tion are shown on the upper right-hand panel of Fig. 1111 
Here the agreement between the BAO measurements using 
the "minimal" and "maximal" angular masks (see Section [5] 
for explanation) is almost ideal. Thus we can conclude that 
our results seem to be robust with respect to the uncertain- 
ties in the survey selection, which is especially true for scales 
with wavenumbers k > 0.03 /iMpc^^. 

The uncertainties in the SDSS photometry might lead 
to the variations in the completeness across the survey area, 
and thus might potentially influence the BAO features. The 
SDSS imaging data is collected along 2.5° wide scan stripes. 
In order to test the effect of varying completeness we have 
created several random catalogues where the completeness 
(and thus the number density of random points) has been 
allowed to fluctuate from stripe to stripe. On the second 
left-hand panel of Fig. [11] we have shown the results of this 
analysis for the case where the completeness was allowed 
to vary uniformly randomly in the range 0.9 — 1.0. The 

Here and in the following we are going to show only FKP 
errors 
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Figure 11. Results of several stability tests. See the main text for a detailed explanation. 
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corresponding right-hand panel shows the case where the 
fluctuations were drawn from the interval 0.75 — 1.0, which 
should be a very conservative choice keeping in mind a su- 
perb photometric quality of the SDSS data. In both cases 
we have plotted ten realizations just to give an idea of the 
possible induced scatter. It is clear from these investigations 
that possible stripe to stripe completeness fluctuations have 
negligible effect on the extracted BAO features. The only 
noticeable variations occur only on the very largest scales, 
k < 0.03/iMpc-\ 

The third topmost panel on the left-hand side shows 
the effect of changing smooth models with respect to which 
the BAO are defined. There we have shown the cases were 
the smooth model is defined by Eqs. (|6]) and (O with the 
"no-wiggle" form of th e linear transfer function given as in 
lEisenstein fc Hiil l| 19981 ). We have plotted the results when 
the smoothing scale a was taken to be a free parameter 
and also when it was fixed to the value of 30 Mpc. The 
third smooth model was defined as a cubic spline fit to 
the measured power spectrum with the spline nodes start- 
ing at wavenumber 0.025 /i Mpc~^ with a uniform spacing 
of O.O5/iMpc-^0 As it was probably expected, with the 
low accuracy of the current BAO measurement the results 
are rather insensitive to how one precisely defines a smooth 
model. 

In our power spectrum calculations we have so far used 
FKP weighting, which is a function of distance only. How- 
ever, with galaxy clusters, in addition to FKP weight, we 
can also apply an extra weight that is proportional to the 
number of member galaxies. This way we give more statisti- 
cal weight to massive systems, and thus one expects to mea- 
sure higher clustering strength. In addition one might expect 
that photometric redshifts are more accurate for the systems 
containing larger number of member galaxies. This is indeed 
what we find: the large-scale clustering amplitude is slightly 
higher and the small-scale damping is mildly weaker com- 
pared to the case with the original weights only. Larger clus- 
tering amplitude along with smaller damping leads to tighter 
errorbars on power spectrum measurement. This can be seen 
on the third right-hand panel of Fig. 1111 Compared to the 
initial weighting the confldence level for the BAO detection 
is increased from 2.2(t to 2.5(t. Apart from slightly different 
confidence for the BAO detection both spectra have oscilla- 
tory features that are in good agreement with each other. 

The lower panels in Fig. [TT] show the effect of changing 
the fiducial cosmological model. For this purpose we have 
generated a MCMC chain with 25, O OP elementiPl for flat 
ACDM models using WMAP^Ifldata (|Noha et alBoogP For 
each of these cosmologies spectra in the parametric form of 
Eqs. ([Sjl and ((Tjl are fitted to the data. Also, as the origi- 
nal power spectrum measurements were done by adopting a 
distance-redshift relation of one particular fiducial cosmol- 
ogy, all of the model spectra are transformed to the "fiducia l 
model framework" in a manner described in iHiitsil l|2006ch . 
The gray "snake" on the lower left-hand panel shows the la 
region filled with the best fitting models. The corresponding 



For this choice of nodes see e.g. iPercival et al. Il2007bl. 
° We have used CosmoMC package l lLewis fc Bridlel2002h avail- 
able at http://cosmologist.info/cosmomc/ 
http://map.gsfc.nasa.gov/ 



right-hand panel shows the histogram of the BAO detection 
confidence level for all of the cosmologies tested. As one 
can see, for a typical flat ACDM model consistent with the 
WMAP5 data BAO is detected at ~ 2cr confldence level. 

There are other possible effects that might influence 
the BAO, e.g. the performance of the cluster finder might 
vary with angular position and distance and thus lead to 
fluctuations in completeness and purity of the final cluster 
sample. However, it is hard to imagine why cluster detector 
should work in a spatially correlated way and imprint a ~ 
100 h^^ Mpc preferable scale to the clustering pattern. Thus 
it is unlikely that the BAO will be noticeably affected by 
those means: even if cluster finding efficiency and fraction 
of "false systems" changes with redshift these changes are 
gradual i.e. can only lead to the changes in the broad-band 
shape of the spectrum. 

In addition, there are always worries related to biasing, 
i.e. how the fluctuating density field as inferred from the 
spatial distribution of the tracer objects is related to the 
underlying matter distribution. Here we would like to argue 
that biasing is probably not a big issue, e.g. if one takes 
N-body simulations of the "concordance" cosmology then 
the typical displacements of the particles are ~ 10 h"^ Mpc. 
These are an order of magnitude smaller than the acoustic 
scale, i.e. in that respect one can look at the structure forma- 
tion as a fairly local process. These small displacements can 
lead only to smooth distortions but cannot produce sharp 
ringing features in the spectrum. 

As an additional test we have investigated the influ- 
ence of the assumed photo-z error model on the detected 
BAO signal. So far in this paper we have used a simple 
model whe re the photo-z errors are independent of redshift, 
whereas in lKoester et al.l (|2007l ) the authors find a mild evo- 
lution where the errors become slightly larger as the dis- 
tance increases. We have used a simple model where the er- 
rors increase linearly with redshift 5z{z) — {z — 0.1)/(0.3 — 
0.1)0.011 + (0.3 - z)/(0.3 - 0.1)0.006, which g ives a good fit 
to the trend observed in lKoester et al.| l|200 7'). It turns out 
that our results in this case are essentially indistinguishable 
from the simpler model assumed earlier. Having a model 
with redshift dependent photo-z errors our model spectra 
are now ca lculated as lig ht-cone averages following the d e- 
scription in iHiitsil (|2006blh (see also l^mamoto et an i|l999^). 
Thus the model spectra are now sums over spectra with dif- 
ferent damping factors. However, for the spatial smoothing 
scales a ~ 30 Mpc the error function in Eq. ((Sjl can be 
simply set to one for wavenumbers k > 0.07/iMpc~^ (with 
an accuracy better than 0.5%), and so our light-cone spec- 
tra are effectively sums over spectra which at most relevant 
k values have approximate power law form, with the only 
difference in their amplitude. Since the sum over power laws 
with varying amplitude and the same spectral index is still a 
power law this should explain why we do not observe any sig- 
nificant deviation from the case when only a single effective 
smoothing scale is used for all the redshifts. Also, the small 
biases of the photo-z errors as measured in iKoester et al-l 
( 2007J turn out to have negligible effect on our final results. 
What concerns the maxBCG clusters then actually we are 
quite lucky as the availability of spectroscopic redshifts for 
many of the galaxies allows one to test the re liability of the 
photo-z errors, and indeed, as seen in Fig. 9 of lKoester et al.l 
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Figure 12. Upper panel: Comparison of the maxBCG effec- 
tive volume with the SPT-like SZ survey covering one octant of 
the sky with the total of ~ 20, 000 galaxy clusters, and with 
the eROSITA-type X-ray survey over 50% of the sky yielding 
~ 100, 000 clusters. Here the two scenarios have been plotted; (i) 
the optimistic case where one has the spectroscopic redshifts for 
all the clusters, (ii) the case with the photo-z errors comparable 
to the accuracy obtained for the maxBCG clusters, i.e. 5z ~ 0.01. 
Lower panel; Contours nP = 3. Below those lines the sampling 
density is high enough for the shot noise not to contribute sig- 
nificantly. The curve for the maxBCG sample is restricted to a 
redshift range z = 0.1 — 0.3. 



(|2007l ) the Gaussian photo-z error model seems to be quite 
adequate. 

In general, for the photo-z errors to have a strong im- 
pact on the BAO demands a nontrivial correlation structure 
over scales of ~ 100^~^Mpc. This might happen at some 
specific redshifts, but it is quite unrealistic for the errors 
to be correlated at those scales throughout the full survey 
volume. 

As a final test we have investigated how the significance 
of the results change for different binnings of the power spec- 
trum. If we use two times smaller wavenumber binning (i.e. 
instead of 16 bins we have now 32 bins) the confidence levels 
for the BAO detection given in Table[2]get replaced by l.Scr, 
2.2(7, and 2. Oct for the FKP, "jackknife", and Monte Carlo 
error models, respectively. The new values for Table |3] are: 
1.6(7, 1.8(7, and l.Scr. 



5 CONCLUSIONS AND DISCUSSION 



troducing a simple model for the nonlinear effects we show 
that "concordance" ACDM model is capable of providing 
very good fit to the data. Moreover, we are able to find 
weak evidence for the acoustic features in the spatial clus- 
tering pattern of the maxBCG sample. If radial smoothing 
scale a, effective bias parameter hes, and nonlinear distor- 
tion parameter q are treated as completely free parameters 
the model with acoustic oscillations is favored by 2.2ct over 
the corresponding "smooth" model without any oscillatory 
behavior. The evidence for BAO weakens somewhat (2.0(7 
in case of Monte Carlo errors), and we are also left with 
some "extra power" on large scales, if we fix the smoothing 
scale (7 to 30 Mpc as suggeste d by the photometric red- 
shift errors 5z = 0.01 estimated in lKoester et al.l (|2007l '). We 
note that these photometric errors actually apply only for 
the very brightest cluster members and in reality the scatter 
for the whole sample might be larger. In addition the clus- 
ter finding algorithm itself might introduce extra smooth- 
ing/overmerging along the line of sight. Due to these uncer- 
tainties we have decided not to carry out any cosmological 
parameter study in this paper, as any untreated systematics 
will be immediately propagated to the parameter estimates, 
leading to biased results. 

We have also carried out power spectrum measure- 
ment by applying additional weights that are proportional 
to the richness of clusters. This way we give more statistical 
weight to systems with a higher mass. As a result we find 
slightly higher large-scale clustering amplitude and mildly 
weaker small-scale damping of the spectrum. The latter ef- 
fect is probably caused by tighter photometric redshift er- 
rors achievable for richer clusters. Larger clustering ampli- 
tude along with smaller damping leads to tighter errorbars 
on power spectrum measurement. Compared to the initial 
weighting the confidence level for the BAO detection is in- 
creased from 2.2(7 to 2.5(7 

E!l 

In our paper we have estimated power spectrum errors 
using three different methods: (i) FKP errors, (ii) "jack- 
knife" errors, obtained by dividing the survey into 5x5x3 
chunks containing equal number of galaxies, (iii) Monte 
Carlo errors found from 1000 mock catalogs. All of these 
errors turn out to agree very well with each other. We have 
tested the stability of our power spectrum measurements by 
relaxing various assumptions made during the main part of 
our analysis. The results of this study, which were presented 
in Section |4l demonstrate the reliability of the original anal- 
ysis. 

The detectability of BAO with ~ 10'' galaxy clus- 
ters is very encouraging result keeping in mind the fu- 
ture cluster surveys, such as the ones based on the mea- 
surement of the thermal Sunyaev-Zeldovich ( SZ) effect 
(jSunvaev fc Zeldovich|[l973 ). e.g. SptE!I, PlanciM or the 
proposed 100, 000-cluster X-ray survey eROSITA^ On the 
upper panel of Fig. [12] we have compared the effective vol- 
ume of the maxBCG cluster sample with the SPT-like SZ 
survey covering one octant of the sky with the total of 



We have calculated the red shift-space power sp ectrum of the 
maxBCG cluster sample (jKoester et al.l |2007| ) - currently 
the largest cluster catalog in existence. After correcting for 
the radial smoothing caused by the photometric redshift er- 
rors, treating the convolution with a survey window, and in- 



This corresponds to the case where the spatial smoothing scale 
due to photometric redshift errors was treated as a free parameter. 

http;//spt. uchicago.edu/ 
^'^ www.rssd.esa.int/Planck/ 
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~ 20, 000 galaxy clusters, and with the eROSITA-type X- 
ray survey over 50% of the sky yielding ~ 100, 000 clus- 
ters. We have plotted two distinct scenarios: (i) the opti- 
mistic case where one has the spectroscopic redshifts for all 
the clusters, (ii) the case with the photo-z errors compa- 
rable to the accuracy obtained for the maxBCG clusters, 
i.e. Sz ~ 0.01. Thus, the eROSITA-like cluster survey with 
spectroscopic cluster redshifts has approximately 50 times 
larger effective volume compared to the maxBCG sample at 
k ~ 0.1 /iMpc^^ implying ~ 7 times smaller errors for P(k). 
If only photometric redshifts with 5z ~ 0.01 are available 
the effective volume increases by a factor of ~ 12 (which 
is also the case for the SPT-like survey with spectroscopic 
redshifts) leading to ~ 3.5 times decrease in power spectrum 
errors. The lower panel of Fig. [T2]shows the contours on the 
wavenumber-redshift plane where nP — 3, which is a typical 
target value of signal to n oise for the fut ure BAO surveys in- 
volving galaxies (see e.g. lAlbrecht et al.ll2006i .l Below those 
lines the sampling density is high enough for the shot noise 
not to contribute significantly. The curve for the maxBCG 
sample is restricted to a redshift range z — 0.1 — 0.3. It 
is evident that except for the lowest redshifts, i.e. z < 0.3, 
those future cluster surveys will have significant levels of dis- 
creteness noise, and thus in terms of the raw performance of 
measuring the BAO only can not compete with a dedicated 
galaxy surveys. 



In reality the main target of these cluster surveys is to 
map the evolution of the cluster number density as a func- 
tion of redshift, since this quantity is very sensitive to the 
amplitude and to the growth rate of the density perturba- 
tions, and as such it provides a powerful tool to constrain 
the properties of DE. The possibility to also find the clus- 
tering power spectrum can be seen as a useful byproduct 
of these surveys that comes essentially "for free" . However, 
both cluster number count study and spatial clustering anal- 
ysis needs as an input estimates for the cluster redshifts, 
with the former probably doing relatively well with rather 
poor redshift estimates (e.g. photometric redshifts). This is 
a rather complicated task since one arguably needs of or- 
der 10 galaxy redshifts per cluster to reliably infer cluster 
redshift. In that respect for the measurement of BAO it is 
certainly less costly to sample the cosmic density field using 
e.g. LRGs or blue emission line galaxies. The latter type of 
objects are targets for th e currently ongoing WiggleZ project 
l|Glazebrook et aLll2007f ) at AAO that attempts to measure 
redshifts for 400, 000 objects over ~ 1000 deg^ and cover 
the redshift range 0.5 < z < 1. Ther e is also a Wfmos in- 
strument (see e.g. iBassett et al.ll2005l l construction planned 
for the Gemini and Subaru observatories that will hopefully 
be completed in 2012. This new multi-object spectrograph 
will be able to measure the spectra of ~ 5000 objects at the 
same time. There are plans to perform a wide field (~ 2000 
deg'^) redshift survey giving spectra for 2 x 10*' galax- 
ies up to redshifts of z ~ 1.3 together with a narrower 
(~ 200 deg'^) and deeper (z ~ 2 — 3) survey with a yield 
of ~ 5 X 10^ galaxies. For the more distant future (around 
2020) one would expect a superb measurement of the mat- 
ter power spectrum using ~ 10® galaxy redshifts obtainable 
via the 21 cm measurements by the Square Kilometre Array 



(SKAj3 (e.g. iBlake et al.ll2004 lAbdaUa fc Rawhndlioosl ). 

However, before Wfmos and Ska become available several 
wide area imaging surveys, such as Pan-STARR^ or Dark 
Energy Survejlff] (Des), will be performed. Although these 
photometric surveys lack accurate redshift information, the 
huge number of objects detectable over large sky areas signif- 
icantly compensate this shortcoming, allowi ng us to obtain 
a highly competitive measurement of BAO ([Blake fc Bridle! 
I2005i ). 

NOTE: 

Several months after the first version of this work was sub- 
mitted a paper appeared ([Estrada et al.|[2009l ) where the au- 
thors perform correlation function analysis using the same 
maxBCG cluster sample. They also find weak evidence for 
the BAO, however their confidence level for the BAO de- 
tection (1.4 — 1.7 (j) is somewhat lower than we typically 
find in our study (1.7 — 2. 2a-). It is imp o rtant to point 
out that similar to our work [Estrada et al.l (|2009l ) consider 
only isotropised quantities. Due to photo-z errors the real 
clustering pattern is significantly anisotropic and it is clear 
that isotropization looses quite some information. Here one 
should realize that in general the information loss for the 
case of the power spectrum somewhat differs from the infor- 
mation loss in the case of the correlation function, or equiv- 
alently, isotropised power spectrum and correlation func- 
tion do not form a Fourier transform pair (only the full 
anisotropic quantities are the Fourier transforms of each 
other). Taking this into account one should not a priori ex- 
pect similar results for the BAO detection. 

It was somewhat surprising to find out that 
lEstrada et all l|2009l l use only five times more random points 
to build their unclustered catalogs. As the cluster sample it- 
self is very sparse this surely leads to a very noisy reference 
level with respect to what the fluctuations are measured, 
leav ing serious doubts on the convergence of their results. 

lEstrada et al] (|2009l ') estimate the correlation function 
covariance matrix in two different ways: (i) "jackknife" 
method, (ii) simple transformation of the power spectrum 
covariance. In their "jackknife" method the subsamples are 
generated by each time randomly removing 1/ 1000th of the 
clusters. This is a valid method only as long as the clus- 
ters are independent. However, due to the large-scale struc- 
ture the clusters are not independent, and thus one should 
use methods which are able to cope better with the situa- 
tion in hand. One method of dealing with dependent data 
is to replace the usual bootstrap /jackknife method wit h 
the block bootstrap /jackknife scheme (see e.g. lLahirill2003l ). 
This is indeed the reason why in large-scale structure stud- 
ies the "jackknife" method is usually implemented by re- 
moving whole contiguous spatial chunks of data, instead of 
just some random set of points. By removing whole spatial 
regions one keeps the correlations inside the regions. If one 
removes random selection of points those correlations are 
lost, which leads to an underestimation of the off-diagonal 
elements in the power spectrum and correlation function co- 
variance matrices. This might explain why their "jackknife" 
error model gives weaker confidence level for the BAO detec- 



http: / /www. skatelescope.org/ 

pan-starrs.ifa.hawaii.edu/ 

https: / /www. darkenergysurvey.org/ 
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tion compared to our case. The other method of estimating 
covariance matrix is actually equivalent to our diagonal FKP 
covariance matrix case. Indeed, their Eq. (17) assumes that 
power spectrum covariance matrix is diagonal. Looking at 
their Table 4 the relevant entry Ax^ ~ 2.4 is quite similar 
to what we find, Ax^ = 2.8 (see Table [5}, and thus is m a 
reasonable agreement with our FKP covariance matrix case. 
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